This document contains additional figures exploring the characteristics of the spawner-recruitment relationships for 267 Pacific salmon stocks located throughout the Northeast Pacific region.
The goal of this document is to examine model performance of various stationary and non-stationary spawner-recruitment model forms for 267 Pacific salmon spawner-recruit time-series spanning five species of anadramous Pacific salmon (Chinook, Chum, Coho, Pink, and Sockeye) spanning from Oregon to Alaska.
## Scale on map varies by more than 10%, scale bar may be inaccurate
#Models
We fit three classes of spawner-recruitment models with 8 models in total, all building off the Ricker stock-recruitment curve:
Dynamic \(\alpha\) (model 3):
\(log(R_t/S_t) = \alpha_t - \beta*S_t + w_t\)
\(\alpha_t = \alpha_{t-1} + v_{t}\)
\(v_t \sim N(0,\sigma_{\alpha})\)
Dynamic \(\beta\) (model 4):
\(log(R_t/S_t) = \alpha - \beta_t*S_t + w_t\)
\(\beta_t = \beta_{t-1} + u_{t}\)
\(u_t \sim N(0,\sigma_{\beta})\)
Dynamic \(\alpha\) and \(\beta\) (model 5):
\(log(R_t/S_t) = \alpha_t - \beta_t*S_t + w_t\)
\(\alpha_t = \alpha_{t-1} + v_{t}\)
\(\beta_t = \beta_{t-1} + u_{t}\)
\(v_t \sim N(0,\sigma_{\alpha})\)
\(u_t \sim N(0,\sigma_{\beta})\)
#Stock Characteristics
Here we include an overview of parameters estimated for each stock fit from the aforementioned models.
Log(\(\alpha\)) estimated from the autocorrelated Ricker function (;model 2) for 267 Pacific salmon stocks. The light histogram is the entire posterior distribution among all stocks, while the dark histogram in the foreground indicates the distribution of median estimates among the 267 stocks:
Summary of the median estimate of log(\(\alpha\)) among all stocks:
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## -0.4601 0.8912 1.2149 1.2683 1.5941 3.0508
Log(\(\alpha\)) estimated from the autocorrelated Ricker function (model 2) broken down by species:
\(Smax\) (1/beta) estimated from the simple Ricker function (model 2) for 267 Pacific salmon stocks:
Smax (1/beta) estimated from the simple Ricker function (model 1) broken down by species:
\(\sigma\) (residual deviance) estimated from the simple Ricker function (model 1) for 267 Pacific salmon stocks:
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## 0.1358 0.5969 0.7940 0.8669 1.0383 3.1441
\(\sigma\) (residual deviance) broken down by species:
\(\rho\), autocorrelation coefficient for residuals from the Ricker function (model 2) for 267 Pacific salmon stocks:
\(\rho\), autocorrelation broken down by species:
#Productivity trends
Time-varying productivity trajectories for each stock, with a 5-year moving average and 95% credible intervals around the average.
### Scaled Umsy